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The dynamics of the n — > 3tt° decay have been studied with the Crystal Ball multiphoton spec- 
trometer and the TAPS calorimeter. Bremsstrahlung photons produced by the 1.5-GeV electron 
beam of the Mainz microtron MAMI-C and tagged by the Glasgow photon spectrometer were used 
for 77-meson production. The analysis of 3 x 10 6 jp — > r\p — * 37r°p — > 67P events yields the value 
a = —0.032 ± 0.003 for the n — * 37r° slope parameter, which agrees with the majority of recent 
experimental results and has the smallest uncertainty. The invariant-mass spectrum was in- 

vestigated for the occurrence of a cusplike structure in the vicinity of the tv + tt~ threshold. The 
observed effect is small and does not affect our measured value for the slope parameter. 

PACS numbers: 14.40.Aq, 12.38.Bx, 13.25.Jx, 25.20.Lj 



I. INTRODUCTION 

The experimental study of the simple and pure strong- 
interaction reaction 



TT TT 



TT TT 



(1) 



is a real challenge as neither a tt° target nor a 7r° beam 
is available. The properties of reaction JT]) can be ex- 
tracted indirectly from complicated processes, for exam- 
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pie, from K + — > TT°TT°e + v e (K^ 4 ), which is the weak de- 
cay of the K + followed by strong final-state interactions 
between the two 7r°s. Major disadvantages of study- 
ing reaction ([1]) in K^ 4 are the small branching ratio 
(2.2 x 10 -5 ) and the complications from the four complex 
form factors for the K^ 4 decay amplitude needed to de- 
scribe the four-particle final state. Another process that 
can be used for the indirect study of reaction (fTJ) is the 
decay 



V 



OTT , 



(2) 



where the 7T°7r final-state interaction can be seen in a 
difference of the r\ — > 37r° decay amplitude from phase 
space. The experimental study of this decay has several 
major advantages: the relatively large branching ratio 
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FIG. 1: Distributions for the phase-space decay of -q — » 37r° (i.e., when a = 0) obtained from a Monte Carlo simulation: (a) 
Dalitz plot, where variable If is the kinetic energy of each of the three pions and (T 1 ) is the mean kinetic energy of the three 
pions (with all energies being calculated in the 77 rest frame); (b) variable z = p 2 /Pmaxi reflecting the density of the Dalitz plot. 



for r\ — > 3ir° (32.5%), a high yield of 77 mesons in many 
production reactions, and very small background from 
other 3ir° contributions, especially in 77 production close 
to the threshold. 

Due to the low energies of the decay 7r°s, 7r°7r° rescat- 
tcring in 77 — > 37r° is expected to be dominated by S and P 
waves. This leads to the parametrization of the 77 — > 3tt° 
decay amplitude as A(j) — > 3ir°) ~ 1 + az p|, where a 
is the quadratic slope parameter that describes the dif- 
ference from phase space. A convenient definition of the 
kinematic variable z is 

3 

z = 6^2(Ei - m v /3) 2 /(m v - 3m ff o) 2 = p 2 /ft 2 na X , (3) 

i=l 

where Ei is the energy of the ith pion in the 77 rest frame, 
and p is the distance from the center of the 77 — > 3ir a 
Dalitz plot. The variable z varies from 0, when all 
three 7T°s have the same energy of m v /3, to 1, when one 
7T° is at rest. A geometrical interpretation of Eq. ([3]) 
gives z = when p = and z = 1 when /? = p max - 
The density of the 77 — > 37r° Dalitz plot is described by 
— + 37r°)| 2 ~ 1 + 2az. The phase-space decay of 
77 — > 37T° (i.e., when a = 0) gives a uniform density of 
the Dalitz plot, as shown in Fig.QJa). The corresponding 
distribution of the variable z is shown in Fig. QJb); it is 
uniform for z from to «0.75 . Experimentally, the slope 
parameter a is usually determined from the deviation of 
the measured z distribution from the corresponding dis- 
tribution obtained by a Monte Carlo simulation in which 
the 77 — > 3ir° decay amplitude is independent of z. 

The 77 — > 37T° decay, which violates G parity, occurs 
mostly because of the u-d quark mass difference. The 
precision measurement of the 77 — -> 3ir° decay width, 
r(?7 — * 3ir°) ~ (rrid — m u ) (1 + 2az), and the param- 
eter a are important tests of Chiral Perturbation The- 
ory (xPTh) 0, In the xPTh momentum expan- 
sion in orders of the xPTh parameter p, the leading 



C(p 2 ) term of the decay amplitude explicitly depends on 
rrid — m u . However, including this term and the second- 
order counter terms, C(p 4 ), is not sufficient to yield 
a decay width that is close to the measured value of 
423 eV PI . The use of dispersion relations [3, H[ > which 
include pion rescattering to all orders, partially improves 
the agreement with the experimental value. For the pa- 
rameter a, the dispersion-relation calculations of Ref. [|[ 
give a negative value in the range —0.007 to —0.0014, de- 
pending on the assumptions made. The results of these 
calculations are outside the value of a = —0.031 ± 0.004 
adopted by the Particle Data Group (PDG) p|. This 
value for a is based on the analysis of 0.9 x 10 6 77 — > 37r° 
decays measured by the Crystal Ball at the AGS la. A 
complete two-loop calculation in standard xPTh [6| re- 
sults in a = 0.013 ± 0.032, the sign of which is even 
opposite to the experimental values. The evaluation of 
the electromagnetic corrections in the 77 — ► 3tt° decay Q 
shows that they are too small to explain the difference 
between the xPTh calculations for a and the experi- 
mental results. Only the use of a chiral unitary ap- 
proach based on the Bethe-Salpeter equation [8| yields 
a = —0.031 ±0.003, which is in very good agreement with 
the PDG value. Several new experiments, which aim to 
remeasure a with better statistics, are still under way. 
So far, the latest preliminary result from the KLOE Col- 
laboration 9], a = -0.027 ± 0.004 sta t ±o.oo6sysu is based 
on poorer statistics (0.65 x 10 6 77 — > 37r° decays) and is 
in agreement with the PDG value within the errors. 

The experimental study of the 77 — > 37r° decay has re- 
cently become of special interest because of new results of 
the NA48/2 Collaboration [l(| that were obtained from 
the analysis of K + — > 7r + 7r 7r decays, where a signif- 
icant cusp effect was observed in the 7r°7r° invariant- 
mass spectrum close to the 7r + 7r~ threshold. The cusp 
occurs because the K + — > 7r + 7r + 7r~ decay contributes 
via the 7r + 7r~ — > 7r°7r° charge exchange reaction to the 
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K + — > 7T + 7r 7r° decay amplitude. Cusps in 7T7r scatter- 
ing were described in the xPTh framework by Meissner 
et al. 11}. The cusp characteristics were used for the 
experimental determination of the rnr scattering length 
combination an — ct2, the xPTh prediction for which is 
0.265 ± 0.004 [l2j. The method for the determination 
of ao — a2 from the analysis of the 7r°7r° invariant-mass 
spectrum from the K + — * 7r + 7r°7r° decays has been pre- 
sented by Cabibbo [l3| . A cusp effect in the 77 — > 3n° 
decay, arising because of the r/ — > 7r + 7r _ 7r° decay contri- 
bution, is expected to be less significant [l4| . This makes 
it less attractive for the experimental extraction of the tttt 
scattering lengths, but neglecting the cusp effect in the 
analysis of the z distribution could result in the wrong 
experimental value for a. In a situation like this, a new, 
high-statistics measurement of the 77 — > 3ir° decays with 
good resolution in the 7r°7r° invariant mass and in the 
variable z is desirable. 

In this paper, we report on a new precision measure- 
ment of the slope parameter a for the r\ — > Sir decay that 
was made by the Crystal Ball Collaboration at MAMI- 
C. These data are also used to look for a cusp structure 
in the 7r°7r° invariant-mass spectrum and for estimating 
how the cusp can affect the result for a. There is also 
an independent analysis of the r\ — > 3ir° data taken by 
the Crystal Ball with a lower beam energy of MAMI- 
ES [Til |l6[ . The result for a reported there is in good 
agreement with the present work. 



II. EXPERIMENTAL SETUP 

The study of the r\ — > 3n° decay was done by measuring 
the process jp — > r\p — > 3ir°p — > 670 with the Crystal Ball 
(CB) multiphoton spectrometer 17J used as the central 
detector and the TAPS calorimeter [l8|, Ell as the for- 
ward detector. The experimental setup was installed in 
the bremsstrahlung photon beam of the Mainz Microtron 
(MAMI) [2^, E2 , with the photon ener gies determined by 
the Glasgow tagging spectrometer [22|, 123, [24[ . 

The Crystal Ball spectrometer was originally buil t by 
SLAC for studies of e + e~ collisions in the J/ip region [25| . 
The recent history of the CB starts at the AGS where the 
first high-statistics measurement of the r\ — > 3ir° slope 
parameter was carried out with the 77-production reaction 
7r~p — ► rjn [5j. After a stint at the AGS, the CB was 
moved in 2002 to the Mainz Microtron for a large variety 
of experiments, including studies of r\ photoproduction 
and rare and special decay modes of -q. 

The Crystal Ball spectrometer is a sphere consisting of 
672 optically insulated Nal(Tl) crystals, shaped as trun- 
cated triangular pyramids, all pointing towards the cen- 
ter of the CB. The crystals are arranged in two hemi- 
spheres that cover 93% of Air steradians. The CB has a 
spherical cavity in the center with radius of 25 cm; it is 
designed to hold a target and inner detectors. There are 
also two tunnels shaped close to a 40° cone; they serve for 
entrance and exit of the beam. Each Nal(Tl) crystal is 



41 cm long, which corresponds to 15.7 radiation lengths. 
As the decay time of Nal(Tl) is about 250 ns, high count 
rates cause pulse pileup and worsen the energy resolu- 
tion. For the runs with the normal count rate, the energy 
resolution for electromagnetic showers in the CB can be 
described as AE/E = 0.020/(^[GeV]) - 36 . Shower di- 
rections are measured with a resolution in 0, the polar 
angle with respect to the beam axis, of ag = 2°-3°, un- 
der the assumption that the photons are produced in the 
center of the CB. The resolution in the azimuthal an- 
gle 4> is cre/smO. More details on the CB spectrometer 
and the physics recently studied with it at the AGS can 
be found in Refs. [M l2Zl El El ■ For experiments at 
MAMI, the CB has been equipped with new electronics, 
providing an individual TDC (time-to-digital converter) 
for each crystal. This gives better suppression of pilcups 
in the analysis of the MAMI data in comparison with 
earlier CB experiments, where one TDC reads out nine 
crystals. 

To cover the downstream beam tunnel of the CB, the 
TAPS calorimeter [U, El was installed 1.5 m down- 
stream of the CB center. TAPS is designed as a versatile 
end-plane hodoscope that can be arranged in different 
formations to optimize the detection of the forward-going 
final-state particles. In this experiment, TAPS was ar- 
ranged in a plane consisting of 384 individual BaF2 coun- 
ters that are hexagonally shaped with an inner diameter 
of 5.9 cm and a length of 25 cm, which corresponds to 12 
radiation lengths. The beam hole of the TAPS calorime- 
ter has a shape of one BaF2 counter removed from the 
hodoscope center. The energy resolution for electromag- 
netic showers in the TAPS calorimeter can be described 
as AE/E = 0.018 + 0.008/(i;[GeV]) - 5 . Because of the 
long distance from the CB, the resolution of TAPS in the 
polar angle 9 was better than 1°. The resolution of TAPS 
in the azimuthal angle (f> is better than 1/R radian, where 
R is the distance in centimeters from the TAPS center to 
the point on the TAPS surface that corresponds to the 
angle. This means that the azimuthal- angle resolution 
of TAPS becomes better than 1° when R > 57 cm. 

The upgraded Mainz Microtron, MAMI-C, is a har- 
monic double-sided electron accelerator with a maximum 
beam energy of 1508 MeV [2l|. The bremsstrahlung pho- 
tons, produced by the electrons in a 10-/xm Cu radiator 
and collimated by a 4-mm-diameter Pb collimator, were 
incident on a 5-cm-long liquid hydrogen (LH2) target lo- 
cated in the center of the CB. The incident photons were 
tagged up to the energy of 1402 MeV using the post- 
bremsstrahlung electrons detected by the Glasgow tag- 
ger [H, El) Ejj- The tagger consists of a momentum- 
dispersing magnetic spectrometer that focuses the elec- 
trons onto the focal plane detector of 353 half-overlapping 
plastic scintillators. The energy resolution of the tagged 
photons is mostly defined by the overlap region of two 
adjacent scintillation counters (a tagger channel) and 
the electron beam energy. For the electron beam of 
1508 MeV, the tagger channel has a width about 2 MeV 
at 1402 MeV (the maximum of the tagging range) and 
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about 4 MeV at 707 MeV (the Ty-production threshold). 
In the analysis, every photon is characterized by the time 
difference (tagging time) between the signal in the corre- 
sponding tagger channel and the experimental trigger. 

The LH2 target is surrounded by a Particle IDentifi- 
cation (PID) detector [3(| that is a 50-cm-long, 12-cm- 
diameter cylinder built of 24 4-mm-thick plastic scintilla- 
tors, which identify charged particles. The PID detector 
was not used in the present analysis. 

The experimental (or DAQ) trigger had two main re- 
quirements. First, the sum of the pulse amplitudes from 
the CB crystals had to exceed the hardware threshold 
that corresponded to an energy deposit of 320 MeV. Sec- 
ond, the number of "hardware" clusters in the CB had to 
be larger than one. In the trigger, a "hardware" cluster 
is a block of 16 adjacent crystals in which at least one 
crystal has an energy deposit larger than 30 MeV. 

A technical paper that will describe more details on 
features, calibrations, and resolutions of our experimen- 
tal setup is in preparation. 

III. DATA ANALYSIS 

The 77 — > 37T° decays were measured by the process 

■yp — > r\p — > 3ir°p — > 672? (4) 

that was extracted by the analysis of events having six 
and seven "software" clusters reconstructed in the CB 
and TAPS together. The six-cluster sample was used 
to search for events of reaction Q in which only six 
photons were detected, and the proton went undetected. 
The seven-cluster sample was used to search for events 
in which all six photons and the proton were detected 
in the CB and TAPS. The cluster algorithm in software 
was optimized for finding a group of adjacent crystals in 
which the energy was deposited by a single-photon elec- 
tromagnetic shower. This algorithm also works well for 
a proton cluster. The software threshold for the cluster 
energy was chosen to be 12 MeV. This value optimizes 
the number of good events reconstructed for the major 
processes: jp — > r\p — > 37r°p, jp — ■> r/p — > 77P, jp — > 7r°p, 
and 7p — » w°ir°p . The hardware read-out thresholds 
for individual crystals were about 1.1 MeV for the CB 
and 3.3 MeV for TAPS. To match the experimental data, 
corresponding software thresholds were introduced in the 
analysis of the Monte Carlo (MC) simulation. 

The kinematic-fitting technique [3l| was used to test 
various reaction hypotheses needed in our analysis. In 
our kinematic fit, the incident photon is parametrized by 
three measured variables: energy and event- vertex co- 
ordinates X and Y in the target. The direction of the 
incident photon is assumed to be parallel to the Z axis. 
The initial values for the mean of X and Y are set to 
zero. The uncertainties ax and cry are taken as the rms 
(root mean square) of the X and Y distribution on the 
target. Their magnitudes, which typically are of a few 
millimeters, are determined by the collimator diameter 



and the 2.5-m distance between the collimator and the 
target. The uncertainty in the energy of the incident 
photon is taken as one-third of the energy width of the 
corresponding tagger channel; this is slightly larger than 
the rms of a uniform distribution. To reproduce the same 
conditions in the analysis of the MC simulation, the sim- 
ulated photon energy is substituted for the energy of the 
corresponding tagger channel. The Z coordinate of the 
vertex is a free variable in the kinematic fit. A photon 
cluster in the fit is parametrized by four measured vari- 
ables. For the CB, they are the cluster energy, the angles 
9 and 4> calculated with respect to the CB center, and the 
effective depth of the shower in the crystals. For TAPS 
instead of the 9 angle, the measured variable is the dis- 
tance R (in centimeters) from the cluster center to the 
Z axis (the same as the TAPS center). In the fitting 
procedure, the 9 angle of a photon is calculated from the 
variable R of the cluster, the distance from the vertex 
to the TAPS surface, and the shower effective depth in 
TAPS. The effective shower depth in the CB and TAPS 
is defined as the distance from the crystal front surface to 
the point where the deposit of the particle energy equals 
half of the cluster energy. The energy dependence of the 
effective depth and its uncertainty for the photon and 
proton clusters were determined from the MC simula- 
tion. 

The resolution of the CB in the cluster angle 9 as a 
function of the cluster energy was determined for photons 
and protons from the difference between initial and re- 
constructed angles in the MC simulation. The resolution 
of the CB in the cluster angle 4> can be obtained from the 
9 resolution by dividing it by sin# (i.e., — ere/ sin0). 
To determine the resolution of TAPS in the distance R 
as a function of the cluster energy, the initial and re- 
constructed values R from the MC simulation were com- 
pared. The resolution of TAPS in the angle <f> can be 
obtained from the R resolution by dividing it by R it- 
self (i.e., 00 = cjr/R). The resolutions so determined 
are used in the kinematic fit as the uncertainties in the 
cluster parameters 9, R, and <f>. 

Besides the individual gain coefficients for each crys- 
tal, which are used to calculate the deposited energy from 
the ADC (analog-to-digital converter) channels, we intro- 
duced an energy-dependent function that provides a cor- 
rection for the cluster energy to get the photon energy. 
The correction is energy dependent bacause less ener- 
getic final-state photons have a larger fraction of the en- 
ergy deposit that is not collected from the crystals due to 
their read-out thresholds. The magnitude of the correc- 
tion was determined from the MC simulation by compar- 
ing the energy of the simulated photons and their recon- 
structed clusters. The advantage of using this function 
is that it improves the invariant-mass resolution and re- 
moves the dependence of the mean invariant-mass value 
on the energy of the incident photon. Without using 
this function, the peaks from ir° and 77 in the invariant- 
mass spectra move to higher masses when increasing the 
incident-photon energy; this is a consequence of more en- 
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ergetic final-state photons. 

The determination of the energy-resolution function 
AE/E for the CB and TAPS was divided into two steps. 
First, we fixed a so-called inherent resolution of the 
calorimeters, which is determined only by the propaga- 
tion of electromagnetic showers in the GEANT simula- 
tion and by the cluster-reconstruction algorithm. The 
inherent energy resolution is better than the experimen- 
tal one, which also includes additional smearing from 
the light collection, the conversion to electronic signals 
in the PMTs (photomultiplier tubes), digitization in the 
ADCs, and pileups. To match the experimental reso- 
lution, which is correlated with the beam intensity and 
depends on the quality of the gain calibrations, the MC 
output for the energy deposited in the calorimeters has to 
be smeared according to an "additional" AE/E function, 
which is different for different experimental conditions. 
For example, to match the experimental energy resolu- 
tion for electromagnetic showers in the CB at a low inten- 
sity of the photon beam (i.e., with few pileups), the en- 
ergy of clusters in the analysis of the MC simulation must 
be smeared according to AE/E = 0.0145/(£[GeV]) - 34 . 
Matching the energy resolution between the experimental 
and MC events can be adjusted via reaching agreement 
in the invariant-mass resolution and in the kinematic-fit 
stretch functions (or pulls) for the energy of the photons 
detected in the CB and TAPS. The superposition of the 
inherent and additional resolution functions is used in 
the kinematic-fit analysis of both the experimental data 
and MC simulation. This resolution function provides 
the uncertainties in the energy of the photon clusters. 
The functions parameterizing the energy resolution for 
electromagnetic showers in the CB and TAPS calorime- 
ter are given in the section describing the experimental 
setup. 

For seven-cluster events, the information from the 
recoil-proton cluster is also used in the kinematic fit. 
However, in contrast to the treatment of a photon clus- 
ter, the energy of the proton cluster is not included in the 
fit, as there is a large uncertainty in the loss of kinetic 
energy of the protons when they travel from the target 
to the calorimeter crystals. This loss occurs in the ma- 
terial located between the target and the crystal surface. 
The calculation of the kinetic energy of the protons from 
the cluster energies is complicated, as the amount of the 
material depends on the production angles in the labo- 
ratory system, and the energy loss in matter depends on 
the kinetic energy of the protons. Moreover, when the 
kinetic energy is above 450 MeV for the CB and above 
370 MeV for TAPS, the recoil protons do not stop in- 
side the crystals. In the analysis of six-cluster events, 
all parameters of the recoil proton (i.e., the kinetic en- 
ergy and the two angles of the momentum vector) are 
free variables of the kinematic fit. The kinematic fit- 
ting includes the four main constraints, which are based 
on the conservation of energy and three-momentum, and 
additional constraints on the invariant masses of certain 
particles in the final state. For our reaction, we can use 



three constraints on the invariant mass of two photons to 
have the 7r°-meson mass and a constraint on the invari- 
ant mass of six photons to have the ?7-meson mass. Then 
the total number of constraints to test the hypothesis of 
reaction ((3]) is eight. The effective number of constraints 
is smaller by the number of free variables of the fit, which 
are the Z coordinate of the vertex and the unknown pro- 
ton parameters (one for seven-cluster events and three 
for six-cluster events). Thus the test of hypothesis Q is 
a 4-C (four effective constraints) fit for six-cluster events 
and a 6-C fit for seven-cluster events. 

Note that testing the jp — > rjp — * 3ir°p — > 6jp hy- 
pothesis means testing this for all possible permutations 
of pairing the six photon clusters to form three 7r°s. For 
six-cluster events, there are 15 such permutations. In the 
case of seven-cluster events, this number is seven times 
larger, as the proton cluster is also involved in the permu- 
tations. The number of permutations can be decreased 
by the separation of the photon clusters from the proton 
one. In our analysis, we used a limit on the 9 angle of 
clusters, which can be only forward ones for the outgoing 
proton, and the information on the time of flight between 
a TAPS cluster and the CB signal with respect to the en- 
ergy of the TAPS cluster. The time-of-flight information 
was also used to remove a small background from the 
six-cluster events, in which one of the clusters was due 
to the proton. The events for which at least one pair- 
ing combination satisfied the hypothesis of reaction ([J| 
at the 2% confidence level, CL, (i.e., with a probability 
greater than 2%) were accepted as r\ — ► 3ir° event can- 
didates. The pairing combination with the largest CL 
was used to reconstruct the kinematics of the reaction. 
A tighter cut on the kinematic-fit CL is unnecessary as 
there is almost no physical background to suppress. A 
looser cut on the CL is not desirable because of inclusion 
of events with poor resolution. 

For the experimental events, in which there are usu- 
ally several tagger hits recorded for one DAQ trigger, 
the 7p — > rjp — > 37r°p — > 6jp hypothesis is tested for ev- 
ery incident photon for which the tagging time is within 
a chosen interval and with the incident-photon energy 
above the reaction threshold of 707 MeV. If an rj — > 3tt° 
event candidate from one trigger passes the 2% CL crite- 
rion for several tagger hits, they are analyzed as separate 
events. The tagging-time spectrum for the experimen- 
tal events selected as r\ — ► 3ir° candidates is shown in 
Fig- [Ha) . The prompt peak from the coincidence of the 
DAQ triggers and the tagged photons sits above a uni- 
form random background. All the experimental spectra 
were produced from the events where the tagging time lay 
between the vertical lines shown in Fig.J^Ja), but are cor- 
rected for the random background by subtraction of the 
spectra produced from the events where the tagging time 
lay outside the prompt region. To decrease the statistical 
uncertainties and fluctuations in the experimental spec- 
tra obtained after the random-background subtraction, 
the width of the random region was taken to be much 
wider than the prompt one. The normalization factor 
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FIG. 2: (a) Tagging-time distribution for the experimental r\ — > Sir events, where the two vertical lines define the prompt 
region, (b) Acceptance for the r\ — > 3tt° events as a function of the incident-photon energy, (c) Measured jp — > r\p excitation 
function compared to the current PWA-fit solution from the SAID [32] data base. 



for the subtraction of the random-background spectra is 
then the ratio of the widths taken for the prompt and 
the random parts of the spectrum. 

Another source of background comes from interactions 
of the incident photons with the target walls. This back- 
ground was investigated by analyzing the data taken with 
an empty target. The size of the so-called "empty-target" 
background depends on the thickness of the target walls 
and on how well the r\ — > 3ir° production on hydrogen can 
be separated from production on heavier nuclei. Since the 
kinematic fit tests the hypothesis of the photon-proton 
interaction, this rejects many r\ — > Sir events produced 
in the target walls. From our analysis, the size of the 
empty-target background in our r\ — ► 3ir° events is 2.2% 
after applying a cut at the 2% CL of the kinematic fit. 
Tightening this cut to the 10% CL decreases this back- 
ground to 1.7%. Another way to suppress the empty- 
target background is to require the detection of the final- 
state proton. In this case, the empty-target background 
in our t] — > 3ir° events is 0.8% after applying a cut on the 
2% CL and 0.65% for the 10% CL The corresponding 
numbers for our experimental 77 — + 37r° events with the 
undetected proton are 7.4% for the 2% CL and 6.0% for 
the 10% CL Because of limited statistics of the empty- 
target data, the remaining empty-target background was 
not subtracted from our experimental 77 — > Sir spectra. 
The results obtained for the slope parameter a by the 
change in the size of the remaining empty-target back- 
ground showed no dependence on it. 

The MC simulation of the 775 — > r\p — > 3ir°p reaction 
was divided into two steps. First, this reaction was sim- 
ulated without dependence of its yield on the incident- 
photon energy. This MC simulation was used to deter- 
mine the 775 — > r/p excitation function. Since we ana- 
lyze the 77? — > r/p — > 3ir°p reaction at a large range of 
incident-photon energies, where the production angular 
distribution changes much, an isotropic distribution was 
used as an input for our MC simulation. The simula- 



tion of the 77 — > 37T° decay was made according to phase 
space (i.e., with the slope parameter a = 0). All MC 
events were propagated through a full GEANT (version 
3.21) simulation of the CB-TAPS detector, folded with 
resolutions of the detectors (such as smearing according 
to the "additional" AE/E function) and conditions of 
the trigger, and analyzed in the same way as the exper- 
imental data. The resulting detector acceptance for the 
775 — > rjp — > 37r°75 events selected by the kinematic fit at 
the 2% CL is shown in Fig. 0Jb); it varies from about 
45% at the 77 threshold to about 25% at an incident- 
photon energy of 1.4 GeV. The experimental yield of the 
IP ~ > VP ~ * 37r°7? events corrected for the acceptance, 
for the 77 — > 37T° branching ratio, for the photon beam 
flux, and for the number of the target protons is com- 
pared in Fig. (2fc) to the current PWA-fit solution taken 
from the SAID [32| data base. The shape of the 77J — > 7775 
excitation function at the production threshold confirms 
the good quality of the energy calibration of the tag- 
ger, which was performed as explained in Ref. [24]. The 
small disagreement that is seen at higher energies can 
be partially due to the difference between the real pro- 
duction angular distribution and the isotropic distribu- 
tion used in our MC simulation. In the second step, the 
775 — > 7775 — > 37r°j? reaction was simulated according to its 
excitation function folded with the bremsstrahlung pho- 
ton distribution. This MC simulation was then used to 
determine the slope parameter a. 

Since the event selection in our analysis is based on 
the confidence level of the kinematic fit, it requires good 
agreement of the \ 2 probability (same as CL) distribu- 
tion for the experimental and MC events. Then a change 
in the cut on the CL value removes the same fraction 
of events from the experimental data and MC simula- 
tion. In Fig. E](a), we compare the CL distributions for 
the experimental (triangles) and MC (circles) events se- 
lected by testing the 77) — > 7773 — > 3ir Q p — > 677) hypothesis. 
These distributions arc in reasonable agreement. Note 



7 



Kin.Fit: y p -> r) p -> 3?t p-»6yp 



s 02 



iSo.15 



0.1 
0.05 




x10° 









[ (a) 


T 


Data 


if 




O 


MC 


s 

" 8 
















i 



0.2 0.4 0.6 0.8 

X 2 probability (C.L.) 



Kin.Fit: 7 p -> 3je p -> 67 p 
x10 6 



S 0-5 



£ 0.4 



0.3 
0.2 
0.1 



(b) 



t Data 

O MC 



5*» 



0.55 0.6 
m(3n°) [GeV/c 2 ] 



Kin.Fit: y p -> r| p -> 3?t p -> 6y p 



in 0.3 



x10° 



c 

UJ 



0.2 



0.1 



(C) 



S 9 



T Data 

O MC 




Kin.Fit: y p -> ti p -> 3?t p -> 6y p 



Constant 


273804.9375 


Mean 


-0.0007 


Sigma 


0.0389 




FIG. 3: (a) The x 2 probability (or CL) distribution for the experimental (triangles) and MC (circles) r\ — > Sir events; (b) 
the 37T° invariant mass for the experimental (triangles) and MC (circles) r\ — > 3tt° events selected at the 2% CL by testing the 
7P — > 3tr°p — > 67P hypothesis; (c) the Z coordinate of the vertex for the experimental (triangles) and MC (circles) 77 — > 37r° 
events; (d) the resolution in the variable z. 



that some increase in the CL distributions for low prob- 
ability is due to events with partially overlapping photon 
showers, or with some leakage of the energy of the show- 
ers from the edge crystals of the CB and TAPS. Since 
the energy-resolution functions for the CB and TAPS 
were determined for "solitary" electromagnetic showers, 
for which all energy is deposited in the calorimeter, the 
errors in the energy and angles are underestimated for 
these "nonideal" clusters. 

The agreement between the experimental data and MC 
simulation for the energy calibration of the calorimeters 
and the invariant-mass resolution can be illustrated by 
a comparison of the 3ir° invariant-mass spectra. These 
spectra obtained from events selected by testing the 
7P — > 3ir°p — > 67P hypothesis are shown in Fig. (3Jb) by 
triangles for the experimental data and by circles for the 
MC simulation. There is good agreement between the ex- 
perimental and MC spectrum for the mean value, which 
is consistent with the 77-meson mass of 547.5 MeV, and 
for the invariant-mass resolution, which has a w 6 MeV. 

Since the Z coordinate of the vertex is a free variable 
in the kinematic fit, the Z distribution must reflect the 
thickness of the LH2 target, which is 5 cm long, and the 



target position. The agreement of these distributions for 
the experimental (triangles) and MC (circles) 77 — > 3tt° 
events can be seen in Fig. (He) . The larger width of the 
Z distribution, compared to the 5-cm thickness of the 
physical target, is due to the resolution of the kinematic 
fit in the Z coordinate. For the jp — > rjp — > 3ir°p — * 6jp 
process, this resolution is about 1.1 cm; it is determined 
from the difference between the initial and reconstructed 
value of Z in the MC simulation. 

The resolution in the variable z, denned in Eq.([3]), can 
be understood from a comparison of z calculated from 
the kinematics of the initially simulated events and z 
calculated from the kinematics reconstructed for these 
events by the kinematic fit. The resolution in z is shown 
in Fig.[3]^d). Based on this resolution and the z-variable 
limits, which are from to 1, we divided our z spectra, 
used for the determination of the slope parameter, into 
20 bins. This binning provides a bin width that is wider 
than one a z — 0.039. The spectrum shown in Fig. [3jd) 
also demonstrates the insignificance of the combinatorial 
background in 77 — > 3ir° — > 67 events (i.e., when the 
kinematic-fit hypothesis with the largest CL chooses a 
false pairing combination of the six photons to the three 
7T°s). The combinatorial background produces accidental 
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FIG. 5: Same as Fig. 0] but for incident-photon energies of 1150 to 1200 MeV. 



reconstructed values for z that are spread widely in the 
dz spectrum, resulting in a washout of the real Dalitz- 
plot slope. For a rough estimate of the combinatorial- 
background contribution, we took the fraction of the MC 
events that have \dz\ > 0.2; this fraction was found to be 
3.5% only. 

To show that our MC simulation satisfactorily repro- 
duces the experimental acceptance, we compare our dif- 
ferential cross sections for the ~/p — > rjp reaction with 
other existing measurements. In Fig. SJa), we show the 
experimental distribution of the rj production angle from 
IP ~ * VP ~ * 37r°p events measured in the center-of- 
mass (cm.) system for incident-photon energies of 800 
to 850 MeV. This distribution is not yet corrected for 
the acceptance. The corresponding acceptance obtained 
from our MC simulation is shown in Fig. SJb) . The com- 
plicated behavior of both the experimental distribution 
and the acceptance is mostly caused by the gap between 
the CB and TAPS calorimeters. The experimental dis- 
tribution corrected for the acceptance, for the rj — * 3tt° 
branching ratio, for the photon beam flux, and for the 



number of the target protons is shown by triangles in 
Fig. 0Jc). This distribution has a smooth shape now, 
which is in good agreement with the jp — > r\p differen- 
tial cross section obtained recently for the same energy 
range by the Crystal Barrel Collaboration at ELSA (CB- 
ELSA) (33|. The CB-ELSA data are shown by circles in 
the same figure. Similarly good agreement with the cor- 
responding CB-ELSA results is observed for all other en- 
ergy intervals. In Fig.[5l we illustrate it just for one more 
energy interval of 1150 to 1200 MeV. The agreement in 
the differential cross sections is shown to demonstrate the 
quality of our r\ — > 37r° analysis. 



IV. DETERMINATION OF THE SLOPE 
PARAMETER a AND ITS UNCERTAINTY 

The experimental z distribution obtained for the r\ — > 
37T° events with prompt tagging times is shown as the 
solid line in Fig. HJa). The events were selected at the 
2% CL from the test of the jp — > r\p — > 37r°p — » 67P 
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hypothesis with the kinematic fit. Both six- and seven- 
cluster events (i.e., without and with the outgoing pro- 
ton detected) were included in this distribution; the frac- 
tion of six-cluster events is only about 20%. The cor- 
responding (unnormalized) z distribution obtained from 
the events with the random tagging times is shown in 
the same figure by the dashed line. The size of this back- 
ground depends on the photon-beam intensity and is less 
than 10% for our data. The experimental z distribution 
after subtraction of the normalized random background 
is shown in Fig. [H^b). The z distribution obtained from 
our MC simulation is shown in Fig. [6jc). This MC sim- 
ulation is based on 4 x 10 7 jp — -> r\p — > 37r°p events with 
r\ decaying to 3ir° according to phase space (i.e., with 
a = 0). The ratio of the experimental z distribution to 
the MC one that was fitted to the function po + 1p\Z is 
shown in Fig.lfJJd). To bring this fit to the required func- 
tion l + 2az, the MC distribution was normalized in such 
a way as to make the fit parameter po equal to 1. Then 
the fit parameter p\ — —0.0322 ± 0.0012 has the meaning 
of the slope parameter a. We take the error of p\ from 
this fit, in which our full set of experimental statistics was 
used, as the statistical uncertainty of the slope parame- 
ter a. To estimate its systematic uncertainty, a variety 
of tests were performed; these are listed in Table |U 
For each test listed in Table U we include information 



on the criteria for event selection, the experimental statis- 
tics, the value for the slope parameter with its statistical 
uncertainty, and the x 2 / ncu? value of the fit. The result of 
the fit shown in Fig. [6](d) is listed in Table [T] as test #1. 
Tests #2 — #4 check the sensitivity of the results to the 
cut on the CL of the kinematic fit. Tightening the cut on 
the CL results in a data set with better resolution and 
less remaining background. The variation in the value 
for a is much smaller than the statistical uncertainties. 

Tests #5 — #7 check the sensitivity of the results to 
a possible background from the 37r° final state that is 
not due to 77 decay. The size of this background can be 
understood from the examination of the 37r° invariant- 
mass spectrum obtained when testing the jp — ► 37r°p 
hypothesis. This spectrum is shown in Fig. [7] for the full 
tagging range 0.7 — 1.4 GeV of incident-photon energies 
and for the ranges 0.7 — 1.1 GeV and 1.1 — 1.4 GeV. An 
examination of the spectrum shows that the size of the 
37T° background comprises a few percent of our full data 
set, and it is negligibly small for events with the incident- 
photon energies below 1.1 GeV. We tested three cuts on 
the incident-photon energy: 0.9, 1.0, and 1.1 GeV. The 
experimental value for a is almost independent of this 
cut. 

Tests #8 and #9 check whether the simulation of the 
threshold on the CB total energy in the DAQ trigger 



10 



Data 


:yp -* 3ji° p, Ey=0.7-1.4 GeV 




30.5 
e 


x10 6 


I Entries 2792804 








a> 




(a) 


5 0.4 








0.3 








0.2 








0.1 








n 




^ i — i — i i 





0.4 0.5 0.6 0.7 0.8 0.9 
m(37i°) [GeV/c 2 ] 




0.4 0.5 0.6 0.7 0.8 0.9 
m(37r°) [GeV/c 2 ] 



0.4 0.5 0.6 0.7 0.8 0.9 
m(37i°) [GeV/c 2 ] 

FIG. 7: Experimental spectra for the 3ir° invariant mass obtained by testing the 7p — > 37r°p — > 6jp hypothesis at the 2% 
CL for (a) the full tagging range of the incident-photon energies, E g — 0.7 — 1.4 GeV, (b) E g = 0.7 — 1.1 GeV, and (c) 



Ea 



1.1 - 1.4 GeV. 



is correct. From the analysis of the sum of the cluster 
energies, the parameters of the trigger were determined 
to be 320 MeV for the threshold itself and a = 20 MeV 
for its uncertainty. To reproduce these trigger conditions 
in the MC analysis, the total energy of clusters has be 
smeared according to a normal distribution having that 
cr, and those events that have a smeared energy less than 
the threshold value must be rejected from the analysis. 
In our tests, we applied software thresholds of 420 and 
470 MeV to both the experimental and MC events. These 
magnitudes were chosen to be considerably larger than 
the threshold of 320 MeV smeared with a = 20 MeV. 
The results obtained for a are in good agreement within 
their statistical uncertainties. 

Tests #10 and #11 check the sensitivity of our results 
to the difference between the isotropic production angu- 
lar distribution used in the MC simulation of the 7p — ► ?/p 
reaction and the real distributions defined by the differen- 
tial cross sections, which depend on the incident-photon 
energy. Examples of the jp — > rjp differential cross sec- 
tions for the two different intervals of the incident-photon 
energy are shown in Figs. [4] and OH In these tests, we de- 
termined the slope parameter from two subsets. The first 
subset included the events with only cos V > and the 
second one with only cos 9 V < 0. Both the results for a 
are in agreement within their statistical uncertainties. 

In tests #12 — #17, we repeated some of the previous 
tests but for the seven-cluster events only. The results 
obtained are in good agreement with each other and with 
the tests that were performed using the sum of the events 
with both the cluster multiplicities. 

Tests #18 — #20 are performed for the six-cluster 
events only. The results for a are slightly smaller here 
in comparison to the corresponding seven-cluster results. 
This could be in part due to a larger fraction of the 
empty-target background in the six-cluster events. The 
rj — > 37T° decay kinematics for this background is smeared 
by the kinematic fit, which assumes the target to be 
a proton. This smearing leads to poorer resolution in 



the variable z and increases the combinatorial back- 
ground, which reduces the real slope of the r\ — > 37r° 
Dalitz plot. A similar smearing occurs for the r\ — > 37r° 
events with accidental incident-photon energies (i.e., our 
random-background events). Test #21, in which a is 
obtained from the r\ — > 37r° events with accidental inci- 
dent photons, illustrates the reduction in the experimen- 
tal slope caused by the smearing effect. Differently from 
the empty-target background, the random background 
was subtracted in all tests of Table HI except test #21. 

Tests #22 — #24 illustrate the stability of the results 
over the period of data taking. Our data set includes 
three periods of data taking from April 2007 to July 2007 
with similar experimental and trigger conditions but dif- 
ferent durations. All three results are in good agreement 
within their statistical uncertainties. 

The uncertainty in measuring the parameter a due to 
our limited resolution in the variable z and the combi- 
natorial background can be estimated by introducing the 
measured a value into the MC simulation and using this 
MC sample instead of the experimental data. To exclude 
statistical fluctuations from this estimate, the same MC 
sample must be used in the ratio of the simulations with 
nonzero and zero a. This can be done when an event 
enters into the z spectrum according to the z value that 
is reconstructed by our program but with a weight of 
1 + 2az, calculated by using z from the simulated kine- 
matics of the event. In this way, the ideal z distributions 
are folded with the experimental acceptance and resolu- 
tion, and both spectra that are used in the ratio have 
correlated statistical fluctuations. These fluctuations are 
canceled in the ratio and do not smear the magnitude of 
the slope. The use of a = —0.032 as an input for this 
estimate resulted in a = —0.0300 ± 0.0007 after applying 
the criteria of test #1 and a = -0.0303 ± 0.0009 after 
applying the criteria of test #4. The difference 0.0017 
between the input a and the value obtained for a using 
the criteria of test #4 is slightly smaller than the value 
obtained using the criteria of test #1. This is expected 
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CL=2%, 7 cl., E 1 < 1.1 GeV 


2.13 x 10 6 


-0.0319 ± 0.0015 
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16 


CL=2%, 7 cl., £ 7 < 1.0 GeV 


1.97 x 10 6 


-0.0319 ± 0.0015 
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17 


CL=2%, 7 cl., £ 7 < 0.9 GeV 


1.62 x 10 6 


-0.0323 ± 0.0017 
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18 


CL=2%, 6 cl. 


0.663 x 10 6 


-0.0292 ± 0.0027 


22.0/18 


19 


CL=10%, 6 cl. 


0.525 x 10 6 


-0.0307 ± 0.0030 
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-0.0301 ± 0.0033 
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because of a tighter quality cut on the CL of the kine- 
matic fit in test #4. In the analysis of the MC simula- 
tion, we can also decrease the combinatorial background 
by applying an additional cut on the difference between 
the initially simulated and reconstructed value of z. The 
elimination of the MC events that have \dz\ > 0.2, ad- 
ditionally to the selection criteria of test #1, results in 
a = —0.0315 ± 0.0007, which is in good agreement with 
a = —0.032 used as an input. These tests illustrate the 
importance of the experimental resolution for the preci- 
sion measurement of the slope parameter a. 

As our final result for the slope parameter a, we use 
the weighted average of all results from Table U except 
tests #18 — #21, which were performed for the six-cluster 
events and the random background. The six-cluster re- 
sults are omitted as they involve larger background and 
much poorer statistics. The weight factor of each result 
was taken as the inverse value of its statistical uncer- 
tainty. This procedure gives —0.0322 for the value of 
a. Note that this value is identical to the result for a 
obtained from test #1, which is based on our full exper- 
imental statistics and could be considered as an alterna- 
tive choice of our main value for a. For the statistical 
uncertainty of a, we take the uncertainty 0.0012 from 
our full-statistics test #1. In the systematic uncertainty 
of a, we include half of the largest variation of the results 



in Table U (except tests #18— #21), which is 0.0013, and 
0.0017 obtained earlier as the difference between the ini- 
tially simulated a — —0.032 and the reconstructed value 
for a. Adding these uncertainties in quadrature gives 
0.0022 for our total systematic uncertainty. 

An uncertainty in the parameter a resulting from a 
possible cusp in the ir°ir° invariant-mass spectrum is 
found to be insignificant. The details on the cusp search 
are given in the next section. 

The value a = -0.0322 ± 0.0012 stat ± 0.0022 syst = 
— 0.0322±0.0025tot, obtained in our measurement for the 
i] — > 37T° slope parameter , is in good agreement with the 
PDG value for a, which is -0.031 ±0.004 Q, but it has a 
smaller uncertainty. That the present PDG value for a is 
identical to the result of the Crystal Ball at the AGS Q 
means that the new Crystal Ball Collaboration at MAMI 
confirms the previous measurement. Finally, taking into 
account the magnitude of our total uncertainty, we prefer 
to give the rounded value 

a = -0.032 ±0.003 tot (5) 

as our final value for the slope parameter. 
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V. SEARCH FOR A CUSP STRUCTURE IN 
THE ttV INVARIANT-MASS SPECTRUM 

The magnitudes of x 2 /ndi for the fits from Table |T] 
hint that the linear hypothesis, 1 + 2az, is not com- 
pletely accurate for describing the slope of the experimen- 
tal r\ — » 37T° Dalitz plot. A possible explanation could be 
the occurrence of a cusp effect in the 7r°7r° invariant-mass 
spectrum from the r\ — > Sir decays. The origin of such a 
cusp was discussed in the Introduction. In this section, 
we check whether a cusplike structure is seen in our data 
and how its presence could affect our result for a. Since 
the effect of the cusp in rj — > 3tt° is expected to be small, 
the most convenient distribution to see it is the ratio of 
the experimental and MC spectrum for the 7r°7r° invari- 
ant mass. This ratio is shown by crosses in Fig. [5{a) for 
the case where the MC simulation of the rj — > 3tt° de- 
cay is made according to phase space (i.e., with a = 0). 
Since the experimental invariant-mass spectrum is 
already distorted from phase space by the nonzero a, 
we also show the ratio in which the MC simulation with 
a = —0.032 is taken instead of the experimental data; it 
is shown by dots in the same figure. This ratio, for the 
most part, is in reasonable agreement with the experi- 
mental one. The largest deviation is seen around a mass 
of 0.28 GeV, which corresponds to the 7r + 7r~ threshold. 
To better separate the effect of the cusp from that of a 
nonzero a value, we replaced the phase-space MC sim- 
ulation in the data-to-MC ratio by the MC simulation 
that includes the slope parameter a = —0.032 for the 
i] — » 37T° decay. The result of this replacement is shown 
by crosses in Pig.[Sfb). 

The magnitude of the cusplike structure, which is seen 
at the level of 1% only, is much smaller than the pre- 
diction made in Ref. [14( (shown by a dashed curve in 
Fig. |HJ). The shape and the sign of the cusplike structure 
also look different. For a more correct comparison, each 
theoretical prediction shown in the figures is folded with 
our experimental resolution. Because of the smallness of 
the effect observed in the 7r°7r° invariant-mass spectrum, 
better statistics is desirable to draw any final conclusion 
on the magnitude and features of the cusp. 

Trying to reach a closer agreement between the data 
and a model prediction, we used formulas from Ref. [34[ . 
The authors of this article conduct a detailed theoret- 
ical study of cusps in Kl — * 37r decays and also indi- 
cate the way to use their approach in the r] —> 3n case 
if one starts from the following parametrization of the 
decay tree amplitudes: A(j) — > 37r°) — uq + u\z and 
A(r/ — > 7r + 7T~7r°) = vq + viy + V2y 2 + V2X 2 with the con- 
ventional rj — > 37r Dalitz plot variables. The cusp shape in 
rj — * 3tt° then depends strongly on the parameters of the 
rj — * 7r + 7r~7r° tree amplitude, which are still not well de- 
termined. To use the parameters available from analyses 
of the density distribution of the rj — > 3n Dalitz plots, we 
assumed that \A(r) -> 37r°) | 2 - 1 + 2(u 1 /u Q )z = 1 + 2az 
and \A(r) -> tt+tt-tt )! 2 - 1 + 2(v 1 /v Q )y + (2v 2 /v + 
vl/vfyy 2 + 2(v 3 /v a )x 2 = l + ay + by 2 + cx 2 . It turned out 



that, after some variation of the rj — > 7r + 7r~7r parameters 
(a, b, and c),by starting from the ones that are available 
in the PDG [l(, the agreement of the cusp prediction with 
the data can be improved. An example of such an im- 
provement is shown by the solid curve in Fig. [5] Further 
improvement requires a simultaneous fit of the r\ — > 'Sir 
and 77 — > 7r + 7T~7r° data, which is not the topic of this 
paper. The current agreement is sufficient to understand 
our systematic uncertainty in the parameter a resulting 
from a possible cusp at the tt + tt~ threshold. In Fig. [9l 
we illustrate the deviation of the slope in the z distribu- 
tion from the phase space for the two cusp shapes that 
were shown in Fig. [5] As shown, if the cusp structure is 
similar to the one that is consistent with our data, then 
the effect of the cusp on the 77 — > 37r° slope is negligible. 
However, the cusp as predicted in Ref. [lj] would cause 
significant distortion away from the linear behavior; this 
is not observed in our experimental data. 



VI. SUMMARY AND CONCLUSIONS 

The dynamics of the rj — > 3ir° decay have been studied 
with the Crystal Ball multiphoton spectrometer and the 
TAPS calorimeter. Bremsstrahlung photons produced 
by the 1.5-GeV electron beam of the Mainz Microtron 
MAMI-C and tagged by the Glasgow photon spectrom- 
eter were used for 77-meson production. The analysis of 
3 x 10 6 7p — > rjp — > 37r°p — > 67^ events yields the value 
a = —0.032 ± 0.003 for the 77 — > 3n° slope parameter, 
which agrees with the majority of recent experimental 
results and has the smallest uncertainty. The agreement 
with a — —0.031 ±0.004 of the measurement made by the 
Crystal Ball at the AGS, where the ir~p — > rjn reaction 
was used, demonstrates the suitability of photoproduc- 
tion for studying subtle effects in rj decays. The 7r°7r° 
invariant-mass spectrum was investigated for the occur- 
rence of a cusplike structure in the vicinity of the 
threshold. The observed effect is small and does not af- 
fect our measured value for the slope parameter. Further 
investigation of the cusp requires better statistics and a 
simultaneous analysis of rj — > Sir and rj — > 7r + 7r~7r data. 



Acknowledgments 

The authors thank J. Gasser and A. Fuhrer for their 
help with the cusp simulation. The support from the ac- 
celerator group of MAMI is much appreciated. We thank 
the undergraduate students of Mount Allison and George 
Washington Universities for their assistance. This work 
was supported by U.S. DOE and NSF (Grant No. PHY 
0652549), EPSRC and STFC of the United Kingdom, the 
Deutsche Forschungsgemeinschaft (SFB 443, SFB/TR 
16), the European Community-Research Infrastructure 
Activity under the FP6 "Structuring the European Re- 
search Area" programme (Hadron Physics, Contract No. 
RII3-CT-2004-506078), DFG-RFBR (Grant No. 09-02- 



13 



MC r) -» 37T, a=0 
O 1 -04 




0.28 0.3 0.320.340.360.38 0.4 

m(7t°7t°) [GeV/c 2 ] 



MC r\ — » 3tc°, g=-0.032 
O 1 -04i 




0.28 0.3 0.320.340.36 0.38 0.4 

m(7c°7c°) [GeV/c 2 ] 



FIG. 8: The ratio of the experimental spectrum for the invariant mass to the corresponding MC spectrum is shown by 

crosses for (a) the phase-space decay of r} — » 3ty° in the MC simulation (i.e., with a = 0) and for (b) the MC simulation of 
the r\ — > 3n° decay with a = —0.032 . The ratio of the MC simulation with a = —0.032 to the MC simulation with a = 
is shown in (a) by dots. The corresponding ratios of the MC spectra with two versions of the cusp prediction are shown by 
dashed (using Ref. [3) and solid (using Ref. [HI) lines in (a) for the case when a = —0.032 and in (b) when a = 0. 
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FIG. 9: Ratio of z distributions (crosses) for the rj — > 3tt° decays simulated with and without a cusp: (a) corresponds to the 
solid curve and (b) to the dashed curve in Fig. |8jb). The solid lines show fits to the function po + 2p\z. 



91330) of Germany and Russia, SNF of Switzerland, and NSERC of Canada. 



[1] C. Amsler el al. (Particle Data Group), Phys. Lett. B 
667, 1 (2008). 

[2] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 539 
(1985). 

[3] J. Kambor et al, Nucl. Phys. B 465, 215 (1996). 
[4] A. Anisovich and H. Leutwyler, Phys. Lett. B 375, 335 
(1996). 

[5] W. B. Tippens et al, Phys. Rev. Lett. 87, 19200 (2001). 
[6] J. Bijnens and K. Ghorbani, J. High Energy Phys. 11, 

030 (2007) rarXiv:0709.0230l [hep-ph]]. 
[7] C. Ditsche, B. K ubis and U. G. M eissner, Eur. Phys. J. 

C 60, 83 (2009) [arXiv:0812.0"344l [hep-phj] . 
[8] B. Borasoy and R. Nissler, Eur. Phys. J. A 26, 383 



(2005). 

[9] F. Ambrosino et al. (KLOE Collaboration), 

larXiv:0707.4l37V l [hep-ex]. 
[10] J. R. Batley et al. (NA48/2 Collaboration), Phys. Lett. 

B 633, 173 (2006). 
[11] U. G. Meissner, G. Muller and S. Steininger, Phys. Lett. 

B 406, 154 (1997) [Erratum-ibid. B 407, 454 (1997)]; 

U. G. Meissner, Nucl. Phys. A 629, 722 (1998). 
[12] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Lett. B 

488, 335 (2000). 
[13] N. Cabibbo, Phys. Rev. Lett. 93, 121801 (2004). 
[14] J. Belina, Diploma thesis, Universitat Bern, 2006, 

|http://www.itp.unibe.ch/diploma_thesis/belina/to talcor.pdf 



14 



[15] M. Unverzagt, Ph.D. thesis, Rheinische Friedrich- 
Wilhelms-Universitat Bonn, Germany, 2007. 

[16] M. Unverzagt et al., Eur. Phys. J. A 39, 169 (2009) 

[17] A. Starostin et al., Phys. Rev. C 64, 055205 (2001). 

[18] R. Novotny, IEEE Trans. Nucl. Sci. 38, 379 (1991). 

[19] A. R. Gabler et al, Nucl. Instrum. Methods A 346, 168 
(1994). 

[20] H. Herminghaus et al, IEEE Trans. Nucl. Sci. 30, 3274 
(1983). 

[21] K.-H. Kaiser et al, Nucl. Instrum. Methods A 593, 159 
(2008). 

[22] I. Anthony et al, Nucl. Instrum. Methods A 301, 230 
(1991). 

[23] S. J. Hall et al, Nucl. Instrum. Methods A 368, 698 
(1996). 

[24] J. C. McGeorge et al, Eur. Phys. J. A 37, 129 (2008). 



[25] M. Oreglia et al., Phys. Rev. D 25, 2259 (1982). 

[26] S. Prakhov et al., Phys. Rev. C 69, 045202 (2004). 

[27] M. E. Sadler et al, Phys. Rev. C 69, 055206 (2004). 

[28] S. Prakhov et al., Phys. Rev. C 69, 042202 (2004). 

[29] S. Prakhov et al, Phys. Rev. C 72, 025201 (2005). 

[30] D. Watts, Proceedings of the 11th International Confer- 
ence on Calorimetry in Particle Physics, Perugia, Italy, 
2004 (World Scientific, 2005) p. 560. 

[31] V. Blobel, Constrained Least Squares and Error Pro- 
pogation, http: / /www .desy.de/~blobel/ 

[32] SAID at http://gwdac.phys.gwu.edu/ 

[33] V. Cred<§ et al. (CB-ELSA Collaboration), Phys. Rev. 
Lett. 94, 012004 (2005). 

[34] M. Bissegger, A. Fuhrer, J. Gasser, B. Kubis and A. 
Rusetsky, Phys. Lett. B 659, 576 (2008). 



